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Abstract — High penetration levels of distributed photovoltaic 
(PV) generation on an electrical distribution circuit present 
several challenges and opportunities for distribution utilities. 
Rapidly varying irradiance conditions may cause voltage sags 
and swells that cannot be compensated by slowly responding 
utility equipment resulting in a degradation of power quality. 
Although not permitted under current standards for intercon- 
nection of distributed generation, fast-reacting, VAR-capable PV 
inverters may provide the necessary reactive power injection 
or consumption to maintain voltage regulation under difficult 
transient conditions. As side benefit, the control of reactive 
power injection at each PV inverter provides an opportunity and 
a new tool for distribution utilities to optimize the performance 
of distribution circuits, e.g. by minimizing thermal losses. We 
discuss and compare via simulation various design options for 
control systems to manage the reactive power generated by 
these inverters. An important design decision that weighs on 
the speed and quality of communication required is whether 
the control should be centralized or distributed (i.e. local). In 
general, we find that local control schemes are capable for 
maintaining voltage within acceptable bounds. We consider the 
benefits of choosing different local variables on which to control 
and how the control system can be continuously tuned between 
robust voltage control, suitable for daytime operation when 
circuit conditions can change rapidly, and loss minimization 
better suited for nighttime operation. 

Key Words: Distributed Generation, Feeder Line, Power Flow, 
Voltage Control, Photovoltaic Power Generation 

I. Introduction 

Displacing fossil-fired generation with renewable genera- 
tion has many desirable outcomes, e.g. reduction in pollution 
and CO2 emissions, and a significant challenge, i.e. reliable 
delivery of electrical power of acceptable quality nearly 
100% of the time OX The mix of renewable generation will 
contain many different resources including wind, concentrat- 
ing solar power, and photovoltaic (PV) at the transmission- 
scale, but with PV as the only presently viable option at the 
distribution scale. The one challenge stated above is actually 
a family of challenges because each of these renewable 
options affects reliability and power quality in different and 
often multiple ways. 

At the transmission scale, renewable generation projects 
are generally large enough to warrant individual transmission 



interconnection studies intended uncover issues that may 
need to be mitigated by the renewable generation owner 
such as installing certain additional equipment or operating 
in certain ways to mitigate the problems. In this case, the cost 
of mitigation is borne by the generator creating the problem. 

At the distribution scale, the size of an individual PV gen- 
erator is so small that the cost of an "interconnection study" 
would be prohibitive. However, when the penetration of PV 
generators on any particular distribution circuit is low, the 
impact is quite small and present utility systems are generally 
unaffected. However, at higher penetrations the net impact 
of many small PV generators may accumulate and affect 
power quality, e.g. slowly responding utility equipment (tap 
changers, switchable capacitors, etc.) not keeping pace with 
cloud-induced rapid variations of PV generation resulting in 
loss of voltage regulation. Fast-response equipment could be 
installed to rectify the problem (e.g. a D-STATCOM Q), 
but the cost is borne by the entire rate base instead of the 
owners of the PV generators who are benefiting from the 
interconnection to the distribution grid. 

A potential solution to the voltage regulation problem 
is to tap into the latent excess PV inverter capacity to 
generate or consume reactive power in an attempt to control 
voltage. Although not permitted by current interconnection 
standards [3], changes to these standards to allow for in- 
jecting or consuming reactive power appear eminent. Under 
this scheme, the burden of providing adequate reactive power 
compensation is again placed upon the generator seeking 
access to the grid. However, many questions still remain 
including: 

• How to dispatch the excess capacity to handle major 
changes in circuit conditions, e.g. rapid change from a 
net real power export to net real power import? 

• How to split the reactive compensation duty equitably 
between the PV generators? 

• Whether the control should be centralized (potentially 
vulnerable), distributed (perhaps more robust), or a 
combination of the two? 

• Whether centralized or distributed, what variables 
should be used as inputs to the control algorithm? 



Despite the challenges related to accommodating PV gener- 
ators, there is also an opportunity for the utility to leverage 
the inverters of these PV generators to enhance its own 
performance such as improving power quality (i.e. voltage 
regulation) and reducing distribution losses via optimal man- 
agement of reactive power flows. However, these should 
be accomplished without placing undue burdens on the PV 
generators by either via excessive dispatch of reactive power 
or by limiting PV generation. 

In the remainder of this paper, we focus on these issues 
surrounding the integration of high penetrations of PV gen- 
eration into distribution circuits. Section [II| discusses these 
issues and their impact on reactive power control systems in 
more detail. In Section [TV) we discuss a few control methods 
that have been proposed and present our own potential control 
method. Section [V] compares these methods via simulations 
of a distribution circuit under widely varying conditions, and 
we draw conclusions and discuss directions for future work 
in Section EH 

We note that there are other approaches to optimizing 
the dispatch of reactive power in distribution circuits for 
the purpose of voltage regulation and loss minimization that 
could be adapted to the present problem including work by 
Baran and Wu H, 0, gl and Baldick and Wu Q and 
also in (SJ and (9). However, these works are somewhat 
specialized to optimal placement, sizing, and/or control a few 
large sources of reactive power where the problem at hand 
includes many small sources. 
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Fig. 1. Diagram and notations for the radial network. Pj and Qj represent 
real and reactive power flowing down the circuit from node j, where Prj and 
Qo represent the power flow from the sub-station, pj and qj correspond to 
the flow of power out of the network at the node j, where the respective 



positive [negative] contributions, p 



and 



and qj 9 ^] represent 



consumption [generation] of power at the node. The node-local control 
parameter q^ can be positive or negative but is bounded in absolute value 
as described in Eq. [4] The apparent power capability of the inverter Sj is 
preset to a value comparable to but larger than maxp^ . 



made either positive or negative. As described in Section IV 
is limited by the apparent power capability of the inverter 



II. Opportunities and Challenges 

Figure [T] introduces the schematic distribution circuit and 
most of the notation we will use in the remainder of the 
manuscript. In our previous work fTOlL ifTTlL we have used the 
LinDistFlow equations ID, (5), (6) to compute voltage and 
power flows. In this work, we switch to an AC solver (T2) to 
compute all distribution circuit quantities. However, we intro- 
duce the problem using the LinDistFlow equations because 
they provide an excellent setting for gaining intuition about 
the competing nature of achieving good power quality, i.e. 
voltage regulation, and reducing distribution circuit losses. 
The LinDistFlow equations for the circuit in Fig. [T] can be 
written 
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where Pj + iQj is the complex power flowing away from 
node j toward node j + 1, Vj is the voltage at node j, Tj +ixj 
is the complex impedance of the link between node j and 
j + and pj + iqj is the complex power extracted at node 
j. Both pj and qj are composed of local consumption minus 
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local generation due to the PV inverter, i.e. pj = py — p 



and qj 
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Pj y) and, q^ are uncontrolled (i.e. driven by consumer load 
and instantaneous PV generation), while the reactive power 
generated by the PV inverter, q^f \ can be adjusted and be 
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Note that reactive power generation is possible only at the 
nodes with PV generation. For the other nodes, we take Sj = 

pf = q f = 0. 

Within the framework of LinDistFlow equations, the rate of 
energy dissipation Cj and the change in voltage AVj between 
nodes j and j + 1 of the distribution circuit are given by 
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= -irjPj+XjQj) 
A. Distribution Loss Reduction vs Power Quality 



(5) 
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Equations ([5]) and ([6]) can be used to discuss many of 
the issues surrounding how to construct a control scheme 
to use the latent reactive power capability of PV inverters to 
maintain power quality and reduce losses. First, Eq. ([5]) shows 
that losses in any circuit segment j are minimized when 
Qj =0. However, to minimize the voltage variation, Eq. ([6]) 
would prefer if Qj = —{rj/xj)Pj in clear competition 
with loss minimization. Therefore, in general, we should not 
expect a control algorithm to simultaneously provide optimal 
voltage regulation and minimize losses. The trade between 
these two desired outcomes must be left up to engineering 
judgement. However, a control scheme should be adaptable 
to easily allow for smooth transitions between emphasis on 
power quality or distribution losses. 



Equation ([6]) also demonstrates the importance of con- 
trolling q^ in a high PV penetration distribution circuit. 
As irradiance conditions change due to cloud passage and 
the pj 9 ^ change rapidly, the segment flows Pj can undergo 
rapid reversals. A distribution circuit that was experiencing 
an acceptable 0.05 p. u. voltage drop without PV generation 
could see rapid switching between the original voltage drop 
and a 0.05 p.u. voltage rise potentially causing voltage 
excursions beyond acceptable bounds. However, if the Qj can 
also be rapidly modified through the q^ 9 \ then the voltage 
variation can be controlled to within acceptable bounds. 

B. Centralized versus Local Control 

Equations ^ and ^ also demonstrate the complexity of 
developing a centralized control scheme versus local control. 
The losses and voltage drop in circuit segment j depend 
upon the flows in segment j, i.e. Pj and Qj. Although not 
currently available to utilities, a centralized controller could 
infer the flows in each segment from smart meter data that 
provides qj and pj for each consumer. With Pj and Qj, a 
centralized controller could determine the dispatch of q^ 
by optimizing an objective function that includes weighted 
measures of losses and voltage deviations. For this type 
of centralized control, the communication requirements and 
additional system vulnerability due to reliance on communi- 
cation may outweigh the potential performance benefits. In 
addition, latency in communication and control may degrade 
performance during rapid changes in cloud cover. 

Local control schemes that act on local variables will not 
suffer from latency and are much less vulnerable as they do 
not depend upon communication for their operation (limited 
communication may be employed by a utility to change 
control algorithms perhaps up to several times during the 
day as overall circuit conditions change fT3ll ). However, truly 
local schemes will only have access to local flows pj and qj 
and, without access to the segment flows Pj and Qj, cannot 
guarantee optimal control. This suggests a local scheme must 
rely upon heuristics to infer enough information about Pj and 
Qj to take appropriate control actions. 

In recent work, ifTOl we have compared centralized and 
local approaches to the control of reactive power. We have 
shown that, for a realistic distribution circuit, a local control 
scheme that simply supplies the local reactive power con- 
sumption (i.e. qj = q^) can achieve almost 80% of savings 
in losses when compared to a centralized control based on 
solving the full optimization problem. Losses were actually 
reduced farther by blending in another heuristic to infer Pj 
and Qj to reduce voltage drops flTTl . The additional heuristic 
works well in reasonably high PV penetration scenarios when 
the circuit is importing or exporting power. However, when 
PV generation and load on the circuit are in relatively close 
balance, the heuristic breaks down and may actually result 
in reduced performance. When a circuit is in balance, the 
Pj randomly change sign from segment to segment, and 
the need to dispatch q^ to regulate voltage is not high. 
Considering the advantages in speed and reliability of local 
versus centralized control and the comparable performance 




Fig. 2. When s is larger than p( g \ the inverter can supply or consume 
reactive power q^ 9 \ The inverter can dispatch q^) quickly (on the cycle- 
to-cycle time scale) providing a mechanism for rapid voltage regulation. As 
the output of the PV panel array p( g ) approaches s, the range of available 
g^) decreases to zero. 

we have simulated in previous work fTOlL fTTlL we only 
consider local control in the remainder of this manuscript. 

C. Equitable Treatment of PV Generators 

Dispatching q^ places additional duty on the inverters 
of individual PV generators which may lead to reduced 
lifetime and increased lifecycle cost. Reference fT3l seeks 
to equitably divide the reactive power duty by setting the 
maximum positive and negative dispatch proportional to 
the capacity of the inverter and PV generator. However, the 
variable that controls q^ between these two extremes is the 
local voltage Vj . Therefore, PV generators that are located on 
a distribution circuit where the voltage is continually above 
or below 1 p.u. will have to endure extra duty compared to 
those located where the voltage is usually close to 1 p.u. 
Retail customers typically have no choice where they are 
located along a circuit. In addition, their location relative to a 
substation may change from day to day depending on how the 
entire distribution system is configured. Therefore, customers 
should not be penalized based on this location. In one the 
alternative control schemes presented in this manuscript, we 
base control of q^ solely on p^ , p^ , and q^ with reactive 

power limits set by the capacity of the inverter so that q^ 
does not depend on location along a circuit. 

In this manuscript, we only discuss the dispatch of reactive 
power and do not consider the calls for limiting PV genera- 
tion. Although we have not encountered situations where this 
control action is required, Ref. fT3ll provides a framework for 
an equitable division of generation reductions. 

III. Modeling Details 

A. Inverter model 

An inverter attached to a PV generator is not an infinite 
source or sink of reactive power. Its instantaneous reactive 
power capability is limited by its fixed apparent power 



capability Sj and the variable real power generation pj . To 
describe this limitation mathematically, we adopt a model 
of PV inverters previously described in ifTOl and iTTH where 
the range o f allowable re active power generation is given by 
\q^\ < \/ s 2 — (p( g }) 2 = q max . This relationship is also 
described by the phasor diagram in Fig. [2] On a clear day 
with the sun angle aligned with the PV array, p^ = pmax 
and the range of available is at a minimum. Although 
Sj relative to p^L could be treated as a free parameter 
subject to optimization, our previous work IfTOl found that 
Sj w 1.1 p^max provides enough freedom in to realize 
the majority of the reduction in distribution losses. Under 
these conditions, \q^\ < 0.45 p$ ax when p^ = p$ ax . 

The choice of Sj « 1.1 Pmax seems reasonable because 
inverters are available in discrete sizes and inverters are likely 
oversized somewhat compared to p^ax 

B. Description of the prototypical distribution circuit 

The configuration of the distribution circuit model we 
consider is similar to one we previously used fTOlL fiTTl . 
Many of the circuit parameters are based on one of the 24 
prototypical distribution circuits described in ifTSll : the nom- 
inal phase-to-neutral voltage Vq is 7.2kV, line impedance is 
(0.33 + 0.3Si)ft/km and constant along the circuit, and the 
distance between neighboring nodes is 0.2 kilometers. The 
circuit consists of 250 nodes, and we study one level of PV 
penetration, i.e. 50% of the nodes include PV generation. 
The capacity of the inverter at each PV-enabled node is set 
to Sj = 2.2 kVA, and the maximum generation capacity is 
set to pmax = 2.0 kW. A uniform level of maximum PV 
power generation assumes identical installations at each PV- 
enabled node (the same pi^lx installed in the same way) and 
spatially uniform solar irradiance. 

We consider two different load/generation cases; under- 
generated and overgenerated. The undergenerated case cor- 
responds to a situation when there is heavy cloud cover 
over the entire circuit and all of the p^ = 0. The load 
at each node is selected from a uniform distribution between 
and 2.5 k W giving an average net real power import per 
node of 1.25 kW. The overgenerated case corresponds to a 
clear sky where all p^ = 2 kW. The load at each node 
is again selected from a uniform distribution, but the limits 
are now and 1 kW . The average generation per node is 
then 1 kW and the average load per node is 0.5 kW giving 
an average net real power export per node of 500 W. The 
reactive power consumed by the loads at each node, q*f \ is 

(c) 

randomly selected from uniform distribution between 0.2p) ' 

(c) 

and 0.3pj corresponding power factors in the range 0.955- 
0.98 which is representative of residential loading (T6). 

The two cases we consider correspond to widely varying 
irradiance and power flow conditions. For a given control 
scheme, the differences between these cases probes the 
robustness of the scheme to rapidly changing irradiance 
conditions. To gauge the sensitivity of the control schemes 
considered in the manuscript to changes in circuit configura- 
tions, we consider many different realizations of a circuit. In 
each realization, the 50% of the nodes that are PV-enabled 



( c ) ( c ) 

are selected randomly and the p K j ' and q K j ' distributions are 
sampled each time. 

C. Power Flow Solution Method 

We use the Matpower package! 12] to solve the AC 
power flow equations. The package implements the Newton- 
Raphson method to solve the equations 

N 

= -(^ c -pf) + ^|^||y fe |(G^cos^A: + ^sin^0) 

k=l 
N 

= -fe c -gf) + ^|^||14|(G^sin^-^cos^) 
k=i 

where dk and Bik are the real and imaginary parts of the 
impedance matrix Y, respectively, and 6^ is the difference 
in voltage phase angle between buses i and k. 

The Matpower package solver does not support situations 
where q\ 9 ^ becomes a function of voltage. For the control 
schemes that require this (described below), we modify the 
implementation of the Newton-Raphson method to take into 
account the variation of q\ 9 ^ with respect to the voltage. 

IV. Control Schemes 

A. Control on Local Voltage Only 

Reference fT3l has proposed a reasonable framework for 
local control of reactive power produced by the inverters 
of PV generators. Although four different modes of control 
are proposed, each consists of a set of piecewise linear 
relationships between q^ and Vj. A simplified version of 
mode 'PV1' from Ref. IU3l is shown in Fig. [3] Although not 
specified in Ref. 11311 , we take q {9) = at Vj = 1 p.u., and 

we take the saturated values of q^ at high and low values 
of Vj to be given by qj 1 ^ defined in the earlier discussion 
of the inverter model. Except for the dynamic definitions of 
qjiax^ fai s j s essentially a proportional control scheme where 

q^ depends linearly on Vj. 

In our attempts to utilize this scheme, the Matpower 
AC solver [12] we employed showed convergence problems 
which we diagnosed as the solution jumping back and forth 
across the points of discontinuous first derivative. To ease 
this difficulty, we smoothed the control function in Fig. [3] 
using a sigmoid function, i.e. 

G(q T - V j9 S) = q- (l - ^ eM J {Vj _ 1)/S] ) • 

(9) 

Here, S is simply a parameter that controls how closely the 
smoothed control function approximates the sharp transitions 
of the original control function. In this work we have taken 
6 = 0.04, and G(qf ax , V, 0.04) is plotted in Fig. [5] for 
comparison to the piecewise continuous control proposed in 
Ref. d. 
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A local control scheme proposed in |[T0l was based on the 
heuristic that losses are minimized when the reactive flows 
Qk are zero, and the were chosen to minimize the net 



Fig. 3. Dark red, piecewise linear curve: A simplified version of the pro- 
posed q^ control function from Ref. fj"3l . Blue smooth curve: Equation fcl 
with S = 0.04 which we use here to improve the convergence properties or 
the AC solver while closely representing the control function in Ref. 1131 . 



B. Control on Local Flows Only 

The control function G(qJ iax ^V J S) describes a form of 
local control that only depends on the voltage at the point 
of inverter connection. Since this voltage is immediately 
available to the inverter, G describes a scheme that would 
be very convenient to implement. However, if the predicted 
replacement of mechanical meters with smart meters occurs, 
information in addition to voltage may be available to control 
an inverter's q^\ We assume a smart meter will be able 
to provide both real and reactive net power flows and that 
these can be communicated to the local PV inverter. The 
inverter will already have measures of its own real and 
reactive power generation. The combination of this data will 
easily provide the inverter with near real-time access to the 

(c) (c) 

three local, uncontrolled power flows, i.e. p K - \ q^, and 

Pj 9 \ It is from these local power flows, as opposed to Vj, 
that we construct an alternative control scheme. We could 
also explicitly include Vj, however as discussed earlier, this 
choice could easily lead to inequities based upon where a PV 
generator is located along a distribution circuit. 

In previous work (TTJ, we have analyzed control schemes 
of the general form 



F k (p 



(9) » (c) « (c) ) 



(10) 



and consistent with constraint ([4]). Here, we summarize 
some of that work. The control scheme is local in that 
q^ depends only on p^ , p^ , q^ . Similar to the voltage 
scheme discussed above, we also assume that the control is 
homogeneous over the line: all inverters are programmed in 
the same way, and explicit dependence on the bus number 
j enters through the inverter's dynamically-determined capa- 
bility q™ ax which in turn depends on Sj and p^ through 
Eq. g). 

It is useful to introduce the following "helper" function, 
Constrj, meant to enforce the constraint ([4}: 
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\q\ < qf a " 
otherwise 



(11) 



reactive power consumption q Q 

Ff^Constr^ 



(ff* at each node: 



(12) 



In this scheme, the inverter supplies the local consumption 
of reactive power up to the limits imposed by its capacity 
Sj and generation pj. This scheme was shown to be very 
effective in reducing the losses. However, as discussed in 
Section [TTJ loss minimization and voltage regulation are 
competing objectives and minimizing losses does not ensure 
voltage regulation. 

In Ref. ifTTIl . the control in Eq. ([12]) was extended to 
consider voltage regulation. Equation ^3]) suggests that, to 
reduce variations in Vj, we should minimize the absolute 
value of the combined power flow VjPj + XjQj. Note that 
for many circuits, the ratio of Tjjxj = a is nearly constant 
for all k and depends only on the configuration and size of 
the conductors used. Thus the absolute value of TjPj +XjQj 
will be exactly zero if for every load node we ensure that 



p 



(c) 



[q^ — Qj 9 ^ — suggesting a control function 

7 aimed at minimizing voltage variations without regard 
for losses: 
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The control in Eq. (12) seeks to minimize losses while 
Eq. ( 13 ) seeks to regulate voltage. A continuous compromise 
between the two objectives in can be achieved via the 
following nonlinear combination 



Fj(K) = Constr fc 



kf; 
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where K is a single parameter controlling the trade off 
between the two objectives in Eq. dTTJ. At K = 1 we recover 
the loss reduction scheme of Eq. (|12|), whereas at K = we 



recover the voltage regulation scheme of Eq. (13). Through 
the parameter K, we now have a simple method to smoothly 
adapt the control scheme, if necessary, as circuit conditions 
change. 

C. Hybrid Control 

We have argued that inclusion of Vj as an input to the 
control method may result in inequitable division of reactive 
power generation duty. However, heuristics used to create the 
control in Eq. ( [T4| ) may, under certain circumstances, fail to 
provide a good estimate of the segment flows Pj and Qj. 
Without knowledge of Vj, the control in Eq. (14) has no 



way of correcting if Vj has moved significantly from 1 p.u. 
To correct this shortcoming, we create a hybrid control that 
combines Eqs. ^ and ( [T4| ). 

The concept behind the hybrid control is similar to that 



used in blend F^ and F^ in Eq. d 14b. We desire that 
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p.u., 



then the control of cn 9 ^ is completely 



governed by Eq. ( 14 ). However, if Vj has fallen significantly 



below 1 p.u., then q 
significantly above 



(9) 



1 p.u., then q)- 



Similarly, if Vj has risen 

_qmax^ s j m pl e 



(9) 



control that achieves this behavior is given by 



H j (K,V j ) = F j (K) + G(q 
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At Vi = 1, G 



and Hj 



F s . For Vj 
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Fj(K) and If, 
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and Hj 



« 1, G -> 
. Finally, if Vj » 1, G -> 
_^max + 2Fj(K)) which 
is still bounded between ±g™ a:c . Here, we have chosen to 
blend the G and F control in one particular way. There are 
clearly other ways to achieve these, but we leave this for 
future study. 

V. Simulations: Results and Discussions 



The control schemes described in Section |IV] are simulated 
on the distribution circuit described in Section IIII-BI The 
node voltages and the distribution circuit losses are calculated 
for both the under and overgenerated cases. For the under- 
generated case, the node voltages are presented in Fig. [4] and 
the losses in Fig. [5] For the overgenerated case, the node 
voltages are presented in Fig. [6] and the losses in Fig. [7] 



A. Base Case-No Control with q 



(9) 



The base case where all q)' 



(9) 



corresponds to the 



situation imposed by the current distributed generation inter- 
connection standards 0. In the undergenerated case where 
Pj and Qj are in the same direction, the voltage deviation 
below 1 p.u. is quite large at about 0.07 p.u. In the 
overgenerated case, Pj and Qj are now in opposite directions 
for the majority of j and, in spite of not taking any actions, 
the maximum voltage rise of about 0.015 p.u. is relatively 
small. In both cases, the maximum deviations take place at 
or very near to the end of the distribution circuit. However, 
during party cloudy daylight hours, the voltage will swing 
0.085 p.u. as the circuit transitions between the under 
and overgenerated cases we consider-uncomfortably close 
to allowable limits. Under higher load or PV generation 
conditions, the voltage swings would easily exceed 0.1 p.u. 
demonstrating the need for control of reactive power in high 
PV penetration scenarios. In the rest of the discussion, we 
use the losses incurred in this base case to normalize the 
losses for the other control schemes. 

B. Control on Local Voltage Only-G(V) 

Controlling the q^ on local voltage via Eq. provides 
excellent voltage regulation with an approximate drop of 
0.027 p.u. in the undergenerated case and a 0.008 p.u. 
rise in the overgenerated case. The total voltage swing on a 
partly cloudy day is reduced to about 0.035 p.u. -a. significant 
improvement over situation when q^ = 0. However, we 
note that the relative losses are increased by about 5% in 
the undergenerated case and by 20% in the overgenerated 
case. The significant increase in the overgenerated case can 
be traced to the rise of Vj over 1 p.u. which forces the 
inverters to consume reactive power increasing the flows 
Qj and the dissipation. The large increase in Qj is in part 



driven by our choice of q 



(9) 




Fig. 4. Undergenerated case: Maximum deviation of Vj from 1 p.u.. 
Black dotted line: No control,g^ = 0; green dashed line: control via local 
voltage, blue solid line: control via local power flows, red dashed line: hybrid 
control. 



significantly from 1 p.u. By reducing these limits, we could 
reduce the dissipation, but voltage regulation will deteriorate 



as the control scheme would begin to resemble q) 



(9) 



0. 



In a comparison of all the schemes (discussed below) and 
in Fig. [8 we show how reducing the q^ limits impacts 
performance. 

C. Control on Local Flows Only-F(K) 

At K = 0, this scheme emphasizes voltage regulation 



j^^rnax 



as Vj deviates 



through Eq. ( [13] ). Therefore, it is not surprising that near K = 
this scheme has similar voltage regulation performance at 
G(V). Near K = 1 where loss reduction is emphasized via 
Eq. (l2\ , F(K) has significantly less dissipation than G(V), 
but the voltage regulation is nearly as poor as q^ = 0. 
Clearly, there is no globally optimum value of K because, 
as we have discussed relative to Eqs. ^ and ([6]), voltage 
regulation and loss reduction are in competition. The choice 
of K is then left up to engineering judgement and Fig. [8] 
(discussed below) provides a useful guide. 

D. Hybrid Control- H(K, V) 

Hybrid control via H(K,V) attempts to contain large 
voltage deviations by smoothly switching from F(K) to 
G(K), i.e. better voltage control, as Vj move significantly 
away from 1 p.u. However, if Vj is close to 1 p.u., H(K, V) 
looks more like F(K) which allows for a greater emphasis 
on loss reduction. For both the over and undergenerated 
cases, the blending of G(K) with F{K) works well for 
voltage regulation with H(K, V) outperforming both G{K) 
and F(K) for K < 1. For losses, the picture is not as clear. 
In the undergenerated case, H(K, V) for K < 1 still results 
in increased losses over the base case and greater losses than 
for G(V). In contrast, H(K,V) for the overgenerated case 
results is significant loss reductions over the base case and 
G(K) around K = 1. The choice of K for this control 
scheme is again left up to engineering judgement. 




Fig. 5. Undergenerated case: Distribution circuit losses normalized by the 



losses when q 



(9) 



0. Lines are the same as in 



it losses 
Fig. g 




Fig. 6. Overgenerated case: Maximum deviation of Vj from 1 p.u.. Lines 
same as in Fig. [4] 
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Fig. 7. Overgenerated case: Distribution circuit losses normalized by the 
losses when qf^ = 0. Lines same as in '. 
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E. Comparison of Control Schemes 

Comparison of the different control schemes is difficult due 
to a lack of a global optimum in both voltage regulation and 
loss minimization-even for a single scheme for a single case 
(i.e. over or undergenerated). In this work, we are interested 
in finding a robust control scheme that can handle the rapid 
variations in power flows as a circuit with a high penetration 
of PV undergoes rapid changes in solar irradiance. Therefore, 
we collapse the over and undergenerated results into a single 
plot by computing the maximum voltage swing experienced 
during the transition from over to undergenerated, i.e. the 
voltages in Fig. [6] minus the voltages in Fig. [4] These values 
make up the vertical axis in Fig. [8] To compare losses, we 
simply average the relative losses from Figs. [5] and [7] and 
these make up the horizontal axis in Fig. [8] Clearly, different 
weightings of the over and undergenerated are possible, as 
are other constructions. However, for this work, we simply 
choose the average. 

In Fig. ^ the points for G(V) and = are two 
ends of what should be a continuous smooth curve because, 
as the limits are decreased, G(V) q^ 9) = 0. The 
point labeled G(V)/2 in Fig. [8] was generated with the limits 
reduced by one half and will also lie on this smooth curve. 
The curve for the control F{K) generally lies higher and to 
the right of the points for G{K) (and its scaled versions) 
showing that F(K) generally gives both poorer voltage 
regulation and higher losses. The F(K) curve does fall lower 
and to the left of the points for G(K) (and its scaled versions) 
for some values of K, but the voltage deviations in those 
cases are already approaching 0.1 p.u. 

The hybrid control H(K, V) and a scaled version labeled 
H(K,V)/2 (where we have reduced the upper and lower 
q^ limits by a factor of 2) generally lie below and to the 
left of G(K) and its scaled versions. The implication is 
that inclusion of local real and reactive power flows into 
a high PV penetration control scheme can lead to better 
performance than simply utilizing the local voltage. The 
curves for H(K, V) and V)/2 are each parameterized 

by K and the family of such curves is parameterized by the 
scaling factor (here we restricted the scaling factors to 1 and 
1/2). Plots, similar to that in Fig. [8] for other distribution 
circuits will clearly aid in the selection of K and the scaling 
factor. 

VI. Conclusions and Path Forward 

In this work, we have developed schemes for controlling 
PV inverter-generated reactive power for high PV penetra- 
tion distribution circuits. In addition, we have developed 
a method for assessing the robustness of these and other 
control schemes during rapid variations in solar irradiance. 
Our metrics of performance included the maximum per unit 
voltage change experienced during the transition from an 
over to undergenerated load condition and the average of the 
dissipation for the two conditions. We have compared the 
control schemes developed in this work to those proposed 
by others fT3l and have reached several conclusions: 
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Fig. 8. Maximum per unit voltage swing experienced on the distribution 
circuit as the loading conditions transition from the over to undergenerated 
case versus the average relative losses in the over and undergenerated cases. 

• The fundamental competition between voltage regula- 
tion and power quality, in general, prohibits control 
schemes from achieving a global optimum, i.e. a mini- 
mum in voltage deviations and circuit dissipation 

• For the cases considered, control schemes that only 
require access to the local variables are sufficient to 
provide adequate voltage regulation. 

• The inclusion of local real and reactive power flows, in 
addition to local voltage, leads to better control system 
performance. 

In this work, we have focused mainly on the rapid tran- 
sitions in loading that a high PV penetration circuit can 
experience during changes in solar irradiance. However, there 
are still open questions related to dispatch of reactive power 
from the PV inverters during other times. For instance, during 
nighttime hours when there is no PV generation and little 
concern about rapid changes in loading, is it equitable to use 
the reactive capability of the PV inverters to improve the cir- 
cuit performance? If so, which control scheme (among those 
considered here or others) provides the best performance? We 
hope that this work spurs others to consider these questions 
in greater detail. 
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